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Abstract 

In the last decade there has been increasing interest in the fields 
of random matrices, interacting particle systems, stochastic growth 
models, and the connections between these areas. For instance, sev- 
eral objects appearing in the limit of large matrices arise also in the 
long time limit for interacting particles and growth models. Examples 
of these are the famous Tracy- Widom distribution functions and the 
Airy2 process. 

The link is however sometimes fragile. For example, the connection 
between the eigenvalues in the Gaussian Orthogonal Ensembles (GOE) 
and growth on a flat substrate is restricted to one-point distribution, 
and the connection breaks down if we consider the joint distributions. 

In this paper we first discuss known relations between random ma- 
trices and the asymmetric exclusion process (and a 2 + 1-dimensional 
extension). Then, we show that the correlation functions of the eigen- 
values of the matrix minors for (3 = 2 Dyson's Brownian motion have, 
when restricted to increasing times and decreasing matrix dimensions, 
the same correlation kernel as in the 2 + 1-dimensional interacting 
particle system under diffusion scaling limit. Finally, we analyze the 
analogous question for a diffusion on (complex) sample covariance ma- 
trices. 
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1 Introduction 



In the seminal paper [I] Baik, Deift, and Johansson prove that the longest 
increasing subsequence of a random permutation has fluctuations governed 
by the (GUE) Tracy- Widom distribution F 2 . This distribution was discov- 
ered by Tracy and Widom as the one describing the fluctuations of the 
largest eigenvalue of random matrices from the Gaussian Unitary Ensem- 
ble (GUE) [38] • 

Soon after, Johansson [26J showed that the same limiting distribution 
occurs in a stochastic growth model, which is equivalent to (a discrete time 
version of) the totally asymmetric simple exclusion process (TASEP) and 
belongs to Kardar-Parisi-Zhang universality class [30] of interacting particle 
systems [33] . This was the beginning of a lot of activities in this field located 
at the intersection between random matrices, stochastic growth models and 
interacting particle systems. For a recent review and a guide to literature 
on the subject, we refer to [2U] or, for a review around TASEP, see |18] . 
Our results are Theorem 11.21 and Theorem 11.31 but before going in further 
explanations, we briefly define the models under consideration. We restrict 
the discussion to the TASEP and Gaussian ensembles of random matrices. 



Gaussian ensembles of random matrices and Dyson's Brownian mo- 
tion (DBM) 

(a) Hermitian matrices. The GUE ensemble of random matrices is defined^ 
as the probability measure on N x N Hermitian matrices H given by 

ih exp {~^ Tl{H2) ) dHi with /3=2 ' (1) 

where dH = Y[^ 1 dH i i Y[i <i< j <N o\Re(Hij)dlm(Hij) and Z N is the normal- 
ization constanto Notice that the measure ([1]) is unitary invariant. 

Dyson [TS] considered a Brownian motion on the space of matrices. More 
precisely, set b id {t) := bj^t) + i&^O), where b} d (t) and t$ d (t), 1 < ij < N, 
are independent standard Brownian motions. The matrix B(t) with entries 
Bi,j(t) := \{pi t j{t) + bj^iit)) is a matrix- valued Brownian motion on Hermitian 
matrices. The stationary matrix-valued Ornstein-Uhlenbeck process defined 

lr rhcre are two other very common standard normalization in random matrix literature, 
see Table 1 in [IS]. 

2 Another way to describe (JTJ) is to take the upper-triangular entries to be independent 
and normal distributed: H ui ~ W(0, N) for i = 1, . . . , N, while Rc(iJ l , J ) - Af(0, N/2) and 
hv.(H itj ) - W(0, N/2) for 1 < i < j < N. 
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by 

dH(t) = -j-H{t)dt + dB(t), with (3 = 2, (2) 

is called = 2 Dyson's Brownian motion and its stationary measure is ([T]). 

(b) Symmetric matrices. The Gaussian Orthogonal Ensemble is a mea- 
sure on N x TV symmetric matrices with probability measure as in ([1]) but 
with /3 = 1 (and, of course, with dH = YLi<i<j<N ^-^ij) ■ Similarly, one 
defines the /? = 1 DBM by (FJ with /3 = 1 (wad ^Jt) = 0). 

Continuous time TASEP 

The continuous time TASEP is a Markov process defined on the space 
Q = {0, 1} Z . For a configuration rj(t) G f2, we say that at position j and 
time t there is a particle if ^(t) = 1, otherwise the position is empty. The 
dynamics is the following: particles jumps to their neighboring right site with 
rate 1, provided the site is empty Let / : fl — > M be a function depending on 
a finite number of rj/s. Then, the backward generator L of TASEP is given 
by 

Lf(v) = 5>a - W d+1 ) - m) (3) 

where if^ +1 is the configuration r) with the occupations at sites j and j + 1 
interchanged. e Lt is the transition probability of the TASEP, see [32],[33] for 
more details on the construction. In the following we will discuss results for 
two specific initial conditions: 

(a) step initial conditions: r]j(0) = 1 for j < and 77^- (0) = for j > 0, 

(b) alternating initial conditions: %(0) = 1 for even j and rjj(0) = for 
odd j. 

In which cases do we have the same limit processes in TASEP and Gaus- 
sian random matrices? The probably most famous result is the convergence 
to the GUE Tracy- Widom distribution F 2 : Let A^™ N the largest eigenvalue 
of GUE N x N matrices. Then [38], 

lim P(A^ x E n < 2N + aN 1 ") = F 2 (s). (4) 

The analogous result for TASEP occurs for step initial conditions^. Let 
x n (t) denote the position at time t of the particle starting from position 
^n(O) = —n. Then, 

lim P(x m (t) < sit/2) 1 / 3 ) = F 2 (s). (5) 

3 A similar result for the partially asymmetric exclusion process has been recently de- 
termined |41j . 
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This connection extends to joint distributions [TTJEZ]- Let ^maxJvW ^ e the 
largest eigenvalue of f3 = 2 Dyson's Brownian motion at time t. Then, 

r ^U 2uN2/3 ) - 2N A , , 

ife = M ^ m 

^ww 3 = Mul 

in the sense of finite-dimensional distributions, where A% is the Airy 2 process 
(firstly obtained in a stochastic growth model by Prahofer and Spohn [36J; 
see also [18] for a definition and properties). 

For GOE matrices, in [39] Tracy and Widom proved that the same rescal- 
ing as in (0J leads to a well-defined limit denoted by Fx and called the GOE 
Tracy- Widom distribution function: 

Mm P(A«° E N < 2N + aN 1 ") = Fx(s). (7) 

JV— >oo ' 

An analogous result holds for TASEP with alternating initial conditions. 
Namely, let x n (t) be the position at time t of the particle starting from 
x n (0) = -2n. Therfl 

lim P(x m (t) < -st^/2) = Fx(s). (8) 

t-»oo 

One might then hope for an extension to the joint distributions of this rela- 
tion. As shown in [T0|l3~7] 

^x,, /4+ „, v ,,W+2 M(t /2)V3 =AM ^ (g) 

in the sense of finite-dimensional distributions, where Ax is the Airy! pro- 
cess (see also [IB] for a definition and properties). However, a convincing 
numerical evidence [7] shows that 

as, aw* — * Mu) - (10) 

The covariance of Ai(u) decays super-exponentially fast in u, while the one 
for the process in l.h.s. of (flOj) only polynomially. 



4 This was proven in [TUM37], but for a related point-to- line last passage percolation 
model, it was obtained by Baik and Rains before [5] , see also [35] for the interpretation as 
growth process on a flat substrate. 
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These results give rise to a number of questions: 

• Is this link between TASEP and Gaussian ensembles of random matri- 
ces just accidential or do they share a common underlying structure? 

• At which point does this variety of connections come to an end? 

The link between GOE and TASEP with alternating initial condition seems 
to be restricted to the static case^|. On the other hand, as we shall discuss 
below, GUE and TASEP with step initial condition have a much stronger 
relation, which can be seen comparing a 2+1 dimensional extension of TASEP 
with the eigenvalues of the GUE minors. However, also this connection is 
only partial: The Markov property at the level of eigenvalues' minor does 
not hold in general as proven in pQ (see also Remark 11.1 in [H]), while it 
holds for the interacting particle system described below. 

2 + 1 dynamics on interlaced particle systems 

An extension of TASEP with step initial condition to a dynamics on a set 
of interlaced particle system has been introduced in [8]. We denote by x™(t) 
the position at time t of the fcth leftmost particle at level m, 1 < k < m < n. 
As initial condition we have x™(0) = k — m — 1 and the configuration space 
of the system with n levels is 

= {x™ e Z | x™ +l <x^< 1 < k < m < n). (11) 

The dynamics is as follows: Each particle x™ has an independent exponential 
clock of rate one, and when the a^-clock rings, the particle attempts to 
jump to the right by one. If at that moment x™ = x™~ 1 — 1, then the 
jump is blocked. If that is not the case, we find the largest c > 1 such that 
x™ = x™^ = ■ ■ ■ = x%£°~l , and all c particles in this string jump to the 
right by one. 

Both the evolution on and its projection onto {x™,m > 1} are 
Markov processes, where the second is nothing else but the TASEP with 
step initial conditions described above. The space-time correlation functions 
for this model are not completely known. However, if we restrict ourselves 
to so-called space-like paths they are determinantal. Introduce the notation 

(ntjti) -< {n 2 M) iff ni < n 2 ,tx >t 2 , and (nx,ti) y£ (n 2 ,t 2 ). (12) 

5 In a related stochastic growth model, the connection extends from the statistics of the 
largest eigenvalue to the one of the the top eigenvalues [17]. However, it is restricted to 
fixed time. 
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We say that (ni,ti) and (n 2 ,t 2 ) are space-like if either (n 2 ,t 2 ) -< (^1,^1) or 
(rii, ti) -< (n 2 , t 2 ). Then, a path is called space-like if any two points on it are 
space-like. The two extreme cases of space-like paths are (1) fixed level n and 
increasing time t and (2) fixed time t and decreasing level n. In [8] it is proven 
that along space-like paths the correlation functions are determinantal. 

Theorem 1.1 (Theorem 1.1 and Proposition 4.2 of [8]). 

For any m = 1, 2, . . . ; pick m (distinct) triples 

Kj = (xj, n^tj) GZxNx R> (13) 

such that 

h < t 2 < ■ ■ ■ < t m , n x > n 2 > ■ ■ ■ > n m . (14) 

Then 

J 3 {For each j = 1, . . . , m there exists a kj, 

1 < kj < nj such that x^.itj) = Xj} = det [/C(xj, ^j)}^ j <m , (15) 

where 

/C(X 1; X 2 ) = -— j ^ dw wXi+ni _ X2 _ n2+l l[(ni,ti)^(n 2 ,t 2 )] 

+ dzi dw 6 [ \ W J -4- , — (16) 

(27ri) 2 J Tl J Tq z w x 1+ n 1+ l gt 2 2(l _ z )n 2 W - Z V ' 

For a set A, Y a is any simple path positively oriented including as only poles 
the elements of the set A. 

Under the diffusion scaling limit 

tf(jTt) - \rt 

one readily obtains that the correlation functions for the A^'s are, along 
space-like paths, still determinantal with kernel 

_ 2 f e (Tl-T2)w 2 -2(£l-&)w 

JC(Ci,ni,r 1 ;C2,n2,T 2 ) = -— / dw; — — l[(m,tiH(n a> t a )] 

Z7ri Jm+e w 

Tll0 2 -25!«) „ .111 1 

dz I dw ^— — (18) 



X fe "(r) := ton ^ (17) 



(27ri) 2 J lzl=£/2 J m+£ e^ 2 ~ 2 ^ z™ w 

where e > is arbitrary. This kernel for T\ = t 2 = 1 appeared first for the 
GUE minors [22] and was shown to occur in TASEP in [TT] . An antisymmet- 
ric version (for general r's) of this kernel was derived in [12] (with a slightly 
different scaling in space) and extends the kernel for the antisymmetric GUE 
minors of I2T1. 
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GUE minors 



The normalization in ([I]) is the best suited to make the comparison clear 
between the scaling limits of large matrices and large time in TASEP. How- 
ever, if we look at matrices of different sizes, it is more natural to drop the 
iV-dependence in the Gaussian term of the GUE measure. Thus, for what 
follows, we consider instead of ((T]) the following probability measure on NxN 
Hermitian matrices: 

J- exp (- Tr(H 2 )) dH. (19) 
Zn 

Denote by A™ the kth smallest eigenvalue of the principal submatrix obtained 
from the first m rows and columns of a GUE matrix. In our context, these 
principal submatrices are usually referred to as minors, and not (as otherwise 
customary) their determinants. The result is well known, see e.g. [6| 1X^122]: 
given the eigenvalues of the NxN matrix, the GUE minors' eigenvalues are 
uniformly distributed on the set 

V {N) = {A™ E R | A™ +1 <\T< Agi 1 , 1 < k < m < N}. (20) 



It is proven in [29J that the correlation functions of these eigenvalues are 
determinantal with correlation kernel 



^ GUE (£i,n i; 6,n 2 ) = -— / dw e ———l 



2m ./ iBJJ w n2- ni 



ni<n2j 



^ / dw-^-— — — , (21) 



(2vri) 2 J VA=£/2 J iR+e e* 2 -^ w - z 1 

for any e > 0. A way of proving is the following. Obviously, changing 
the condition A^ +1 < A™ into A™ +1 < A™ does not change the system, 
since we cut out null sets. Then, using Sasamoto's trick originally employed 
for TASEP [37], one can replace the interlacing condition by a product of 
determinants, 

JV-l 

ndet[0(Ar,A J m+1 )] 1 <„< m+1 , (22) 

m=l 

where A™ +1 = virt are virtual variables, (f>(x,y) = l[ x < y ], 0(virt,y) = 1. 
Thus, the measure on T>^ N ' becomes 

(JV-l \ N 

H det[0(Ar,Af +1 )] 1 < M < m+1 A(X N )l[e-^ 2 dA, (23) 
m=l / i=l 

where dA = Yli<k<n<N and A is the Vandermonde determinant. Finally 
one simply applies Lemma 3.4 of [10]. A further approach is presented in |22j . 
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Results 

Evolution of GUE minors 

Does there exist an extension of such a result for the minors of (3 = 2 Dyson's 
Brownian motion? There are two aspects to be considered. The first is to 
determine whether the evolution of the minors' eigenvalues can be described 
by a Markov process. It is known that it is not the case if one takes at least 
three consecutive minors pQ . However, along space-like paths the evolution is 
indeed Markovian (see Section H]). The second issue concerns the correlation 
functions and if they have any similarities with the ones for the 2 + 1 particle 
system defined above. As we shall prove, the answer is affirmative if we 
restrict ourselves to space-like paths. This is the content of Theorem 11.21 
below. 

To make the connection more straightforward, we replace the Ornstein- 
Uhlenbeck processes by Brownian motions starting from 0. Note that the 
two models are the same after an appropriate change of scale in space-time. 
While preparing this manuscript the result analogue to Theorem 11.21 below 
for the Ornstein-Uhlenbeck case was obtained by Adler, Nordenstam and 
van Moerbeke [2]. 

Let H(t) be an N x N Hermitian matrix defined by 



where b^j{t) and bij(t) are independent standard Brownian motions. The 
measure on the N x N matrix at time t is then given by 



For n G {1, . . . , N} we denote by H(n, t) the n x n minor of Hit), which 
is obtained by keeping the first n rows and columns of H(t). Denote by 
^i(t) < < • • • < A™(£) the eigenvalues of H(n,t). Then, at any time 

t, the interlacing property ( 12 Up holds. Moreover, along space-like paths the 
eigenvalues' process is Markovian with correlation functions given as follows. 

Theorem 1.2. For any m = 1,2, ... , pick m (distinct) triples 




(24) 




(25) 



>o 



(26) 



such that 



ti<t 2 < 



■ < t 



n\ > n2 > ■ ■ ■ > n, 



(27) 



S 



Then, the m-point correlation function of the eigenvalues' point process is 
given by 

p^(xi, . . . , x m ) = det [/C GUE (x,, x.OU^, (28) 

w/iere 

2 /" e (ti-*2)«' 2 -2(a;i-X2)io 

/C GUE (x i; x 2 ) = - — J l[KAH(n 2 ,t 2 )] 



2 T f e w 2 ti-2xiw j 
+ 7 rrr <P dz dw — ; (29) 

(27Ti) 2 J[,| =£/2 y iR+e e ^" 2 ^ w-zz^ 



where e > 0. 



Evolution of Wishart minors 

The appearence of determinantal correlation functions along space-like paths 
is not only limited to Brownian motion on GUE matrices, but they also occur 
in other Hermitian matrix models, namely the Laguerre Unitary Ensemble. 
We show that the evolution of Wishart matrices [40J along space-like paths 
is determinantal and determine the space-time correlation kernel. 
Let A(n, t) be a p x n complex valued matrix defined by 

A itj (n,t) = -/=(M<) 1 < <P, 1 < i < « (30) 

where the b^/s and bij's are independent standard Brownian motions. Then, 
we define the (complex) n x n sample covariance matrix (or Wishart matrix) 
by H(n,t) = A(n,i)*A(n,i), which is usually referred to as the Laguerre 
process. As before, denote by A£(i) the fcth smallest eigenvalue of H(n,t). 

Theorem 1.3. For any m = 1, 2, . . . , pick m (distinct) triples 

Kj = (xj,nj,tj) e R x {1, . . . ,p} x R> (31) 

such that 

h <t 2 < • ■ ■ < t m , n l >n 2 >---> n m . (32) 

Then, the m-point correlation function of the eigenvalues' point process is 
given by 

p™{x x , . . . , x m ) = det [/C LUE (x 4 , M j )) 1 ^ m , (33) 

where 



1 /• „xi/(z-ii) _ + \p-l-m 



K (X i; X 2 ) = -— ^ dz (2 _ t2)p+1 _ n2 l[(n ll t 1 )^(n 2 ,t 2 )] 

-i r r e *2/(z-ti) t z _ t w p i 

dz ® dw 77 rr^ 7—— . (34) 



(27ri) 2 / ro J rz ( e xa/(w-t 3 ) ( w _ £ 2 ) P +i-n 2 z p ^ _ z ' 
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2 Proof of Theorem EH 

The issue of the Markov property is discussed in Section|4]below and therefore 
we assume it to hold in this section. For < ti < t 2 , the joint distribution 
of Hi = H(n,ti) and H 2 = H(n,t 2 ) is given by 

const x exp exp f_ M&Zm \ <iHl <m . (35) 



The measure on eigenvalues is obtained using Eynard-Mehta formula 
for coupled random matrices, which on its turn is based on the Harish- 
Chandra/Itzykson-Zuber formula [231 El] (see Appendix [D]). It results in 
the following formula. 

Lemma 2.1. Let n be fixed. Denote by A£(i), 1 < k < n, the eigenvalues of 
H(n,t). Their joint distribution at < t\ < t 2 is given by 

const x A(A n (ti)) det ( e -W(*0-*?fcO) a /(fa-*iA A(A n (t 2 )) 

V / l<i,j<n 

N 

x J] e -W(*0) 3 Ai dA »( fl ) dA"(t 2 ), (36) 



i=l 



with A the Vandermonde determinant and \ n (t) = (A™(t), . . . , A™(t)). 

The second formula concerns the joint distribution of the eigenvalues 
at two different levels. This result is a special case of the formula f!23p 
discussed above. (It is enough to reintegrate out the lower levels, which gives 
a Vandermonde determinant). 

Lemma 2.2. Let t be fixed. Denote by AjJ(t), 1 < k < n, the eigenvalues of 
H(n,t). Their joint distribution at levels n and n + 1 is given by 



const x A(A"(t))det[0(AI l (t),A J n+1 (t))] 1 <„< n+1 A(A" +1 (t)) 

n+1 

x J] e-C^W/t dA^(t) dA" +1 (t), (37) 



i=l 



where A^ +1 = virt are virtual variables, 4>(x,y) = l[ x <y], 0(virt, y) = 1 (and 
A the Vandermonde determinant) . 
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The eigenvalues' process is a Markov process (see Section H] for details) 
for both fixed matrix dimension n and increasing time t, as well as for fixed 
time t and decreasing matrix dimension n. The combination of the formulas 
in Lemma 12.11 and Lemma 12.21 leads to Proposition 12.31 

Proposition 2.3. Let N± > • • • > N m = 1 be integers and < t\ < ■ ■ ■ < t m 

be reals. We denote by A™(t) < • • • < A™(£) the eigenvalues of H(n,t) and set 
N = Ni, N m+ i = 0. Then the joint density of 

{Xk(tj) : 1 < j < m, Nj < n < Nj-i, 1 < k < n} (38) 

is given by 
const xdetK^^fCii))]^^ 

det [T tj+1 , tj (Af J (t j+1 ), Af (tj))] 1<k4<N4 



><n 



iV 

X 



n det[0(Ar 1 (t, +1 ),A^(t J+1 ))] 1 < M < n 

n=AT i+1 + l 



(39) 



where 



<i>(?,y) = l[ t < v ], (j){xl ,2/) = l, 

7t, s (x, y = exp l [t > s] , 

- s) V t-s / (40) 

/or fc = 0, . . . , JVi — 1. Here is the standard Hermite polynomial of degree 
k (see Appendixllft for details). 

We could have chosen any polynomials of degree k multiplied by the 
Gaussian weight without changing the probability measure ( 1391) since the 
modifications would just affect the normalization constant. However, this 
choice allows a huge simplification of the computations, because of the prop- 
erties of Lemma 12.41 below. 

To determine the kernel, we first slightly rewrite (13"9"|) . For 1 < n < N\ 
let c(n) = #{z : Ni = n}, and we denote the consecutive times for such a 
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level by t" < ■ ■ ■ < t", n . Then, the measure (I3"9j) can be rewritten as 



c(n) 
Ni / 

const x J] det[0(Ar 1 (tr 1 ),AF(C (n) ))] 



l<fe,^<n 



n=2 

c(n) 

X 

o=2 



Ildetfe^^A^C), A?(Ci))]^ B ) det[<i (Af (ff 1 ))]^^. 



(41) 



It is known that a measure of this form has determinantal correlations and 
the correlation kernel is computed by means of Theorem 4.2 of [9], which we 
report in Appendix [A] for the reader. 
For any given fcGZwe set 

*r<*> = lj^"- 2W "^ k > «>0- (42) 

For n — N\, t — t\ and fc = 0, ...,N\ — 1, this function is the one in 
the measure f[3"9l . which is obtained from the first representation of Hermite 
polynomials in fill 01) . 

Lemma 2.4. It holds, for < r < s < t and k > 1, 

(i) * K'- k = K-d k , 
(a) Tt, s * K'\ = K't k , 

(Hi) <f)*Tt >s = Tt, s * 4>, 

(iv) % >s * T s>u = T t>u . 

Proof. For the first relation, we use Re(w) = e > so that we can exchange 
the two integrals, 

l( ,w 2 ~2wy/t 1 / 2 w k 



poo 2,^+1 ■y r 

= J m tytmr^L^ 

P POO 

/ dwe w2 w k / dye- 2wy/tl/2 (43) 

JiR+e Jx 



f (fc+i)/2 2vri 
2 k 1 



t k l 2 2vri 



./iR+e 
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For the second identity, we first do the change of variable w = z(s/t) l l 2 
in the integral representation (|42|) of and then perform a Gaussian in- 
tegration: 

t (K+l)/A 27T1 J iR+£ J R / n (t - s) 



2 k+l 1 



(44) 

The third relation is also easy to verify. Indeed, 

i>*T t , s )(x,z) = / dy(f)(x,y)T ttS (y, z) = / dy T t>s (y + x , z) 



dy%, s (x,z-y)(j)(z-y,z) = / dyT t , s (x, y)<j)(y, z) = (T t , 8 * z). 

(45) 

The last relation is the standard heat kernel semigroup identity. □ 

By Theorem IA.1I and Remark \A.2\ there is a simple way of getting the 
kernel if the matrix M with 

m m = (</> * r fc * ■ • • * <t> * r^ 1 * ^i 1 )^- 1 ), (46) 

is upper triangular, where T n := %« t t n - The identities in Lemma [2.41 give, 
for k > e, 



Vo' for£ = fc (47) 

because the last expression is (after a rescaling in x) proportional to the 
orthogonal relation fllQ9p for n = and m = k — i. 

Next we need to determine the polynomials $"'*(a;), I — 0, . . . , n— 1, which 
are biorthogonal to the functions ^^(x), fc = 0, . . . , n — 1, i.e., polynomials 
satisfying 

f dx y n k '\x)<S> n £ '\x) = S k/ , 1 < k, I < n - 1. (48) 
Lemma 2.5. T7ie functions 

= iVwf^) = # dz 6 ^ (49) 



?! 2* ^ 2* 2yri 

satisfy the relation 



7r 
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Proof. One does the change of variable x i— > x\/t and then uses the orthog- 
onal relation (1X05]) . □ 

Let us compute the last term in fllOSp . To simplify the notations, we set 
t\ = t™* and t 2 = t™ 2 . Fhst , we do the changes of variables w = and 
z = ^T 2 l in (^D and (in]). We obtain 



n 2 n2 nni OP r p w 2 t\—'lwx\ ~.n\~k 

E KU^TM) = E ^ p-p « £ n V „_ 2& , ,„, +1 _ t 

(50) 

Now, we take the integral over 5 to satisfy |z| < |u>|, say \z\ = e/2. This 
allows us to take the sum inside and extend it to +00 (because for k > n 2 
the pole at zero for z vanishes). The sum over k gives 



~ z k-i 



(51) 



w — z 

k>l 



so that we obtain 

"2 



k=i 

r,2 



= — 7 — / dw ffi dz — ^- — — . (52) 

2-2 (2vri) 2 7 iR+e Jj a | =£/2 e - 2 *2-2-2 z m w - z 

The last term we have to compute is ( - ta i' ta 2^ see Theorem 11.11 To 
simplify the notations, we set <p <yta i' ta, 2\x,y) = <p <yni,tl ' ,fl2,t2 \x,y). We have 



(ni,ti;ri2,t2) 



* (n2 - ni) *T t2 , tl , if KtiH (n 2 ,t 2 ), 
0, otherwise. 



(53) 



It is easy to verify that 0(x, y) has the integral representation 



and similarly, 



2 f e 2w(y-x) 
27T1 l im , „ 2w 



9 /" p 2w(y-a;) 

(f)* n (x,y) = / dw- — , e>0. (55) 



Then, for (ni,ii) -< (^2,^2), a Gaussian integration gives us 

oni 9 /" p w 2 (t\-t2)-2w(xi-X2) 

^ 1 2; 2«2 2vri J iR+e w n *-^ 1 ; 
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Equations f[52~j) and (|5l)|) yield a kernel which is, up to the conjugation 
factocl 2" 1- " 2 , the same as (T2"9"j) . Thus the proof of Theorem II .21 is completed. 

3 Proof of Theorem 11.3 



As for the GUE case, the issue of the Markov property is discussed in Sec- 
tion H] below. For < ti < t 2 , the joint distribution of A\ = A{n,t\) and 
A 2 = A(n, t 2 ) is given by 

/ Tr(AMi)\ / Tr((A* 2 -A* 1 )(A 2 -A 1 ))\ t 1A . , 
constxexpf y 11 j exp ( vv 2 ^ ^ 2 ^JdAidAs. (57) 

The measure on eigenvalues is obtained (as in the Ornstein-Uhlenbeck case 
studied in [10]) by the Harish-Chandra/Itzykson-Zuber formula for rectan- 
gular matrices j25jH2] (see Appendix |D]) . It results in the following formula. 

Lemma 3.1. Let n be fixed. Denote by A£(£) , 1 < k < n < p, the eigenvalues 
of H(n,t) = A(n,t)*A(n,t). Their joint distribution at < t\ < t 2 is given 
by 



const 
x det 



/ VWAp\ / A?(fa) x (p-n)/2 {xnti)+X7{t2))/{t2 _ ti) 

L v- n \ t 2 -h j Uf(ti) J e 

\ / v 7 J l<i,j<n 

n 

A(A"(t 1 ))A(A"(t 2 )) n(Ai(^i))^ n e- A?( * l)Al dA?(tO dA?(* 2 ), (58) 
i=i 



where I m is the modified Bessel function of order m, see ftllty. 



The second formula concerns the joint distributions of the eigenvalues at 
two different levels. This is studied in [22] with the following result. 

Lemma 3.2. Lett be fixed. Denote by A^(t), 1 < k < n < p, the eigenvalues 
of H(n,t). Their joint distribution at levels n and n + 1 is given by 

const x A(A^(t))det[0(Ar(t),A J "(t))] 1 ^.< ri+1 A(A" +1 (t)) 

n+l 

x n{K + \t)) P ~ {n+1)e ~ K+1{t)/t dK(t)d\? + \t), (59) 
t=i 

where A™ +1 = virt are virtual variables, 4>(x,y) = t\ x >y] and 0(virt,?/) = 1. 



6 A determinantal point process is defined by its correlation kernel, which is defined up 
to conjugations. 
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Putting together the formulas in lemmata 13.11 and 13.21 leads to the next 
proposition. 

Proposition 3.3. Let p > N% > ■ • • > N m = 1 be integers and 
< ti < ■ ■ ■ < t m be real numbers. We denote by \™(t) < ■ ■ ■ < A"(i) the 
eigenvalues of H(n,t) and set N = iV 1; iV m+1 = 0. Then the joint density 
of 

{Xk{tj) : 1 < j <m,Nj <n< N^, l<k<n} 



(60) 



is given by 
const xdetKy(A* (h))] 

m— 1 - 

n det[7j;3(Af ( tj+1 ),xf j 



X 



<kl<rii 



i=i 



AT,- 



X 



det[0(Ar 1 (t J+1 ),A^(t J+1 ))] 



Kfc,^<ri 



(61) 



where 



(x,y) = l [x > y] and (f)(\™ +1 ,y) 



V! s (x,y) 



n/2 



2^/xy 



1 



exp 



t — s / t — s 

, p-ATi 



X + ?/ 

't-s 



l[a;,j/>0]l[s<t]) 



X 



(p - N x + k)\t\ +1 V*i 



exp 



X 

— 1 llx>0]j 



(62) 

for k — 0, . . . , N\ — 1. i/ere L£ are i/ie generalized Laguerre polynomials of 
order n and degree k, see AppendixXCk 



Comparing the mathematical structure of fl39|) and fjSTj) . we see that the 
only difference is that the transition kernel for time depends also on the level. 
However, this does not pose any problem, see Remark IA.31 

For k G Z and x G IR+ we set 



j-(*+i) 
2?ri 



dz 



a x(*-l)/t 



r 



n+fc+l 



(63) 



For n — p — Ni, t — ti and fc = 0, . . . , N\ — 1 the above defined function coin- 
cides with fl6"2~|) . Moreover, the prefactors are chosen such that the following 
nice recursion relations hold. 
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Lemma 3.4. It holds, for t>s>r>0,n<p, and k > 1 

(a) vr * KT = K-V, 

p—n ,y-p— (n— 1) , 



(Hi) (p*l7 s " = 7? s 



To prove this lemma, we first obtain a different integral representation 
for (jH3D- Namely, after the change of variable z = z/(z — t) we get 

mf( x ) = zl I & ( z -i)"' 1 e *az-t) (64) 

k v ' 2vri/ ro z n+k+l y ' 

Proof of Lemma \3.4\ Using the representation (|64|) . we have 



r 

— / dz-^=_2! (e*/(*-*) - 1) = ^rf'Va;) 



(65) 



because for k > the term independent of x has residue equal to zero. 

Using the integral representation (j!17j) of the modified Bessel function I n 
in flBTj) . we get (for x, y > 0, t > s > 0) 

^y) = o-:/, 3 <P 3nT«P ~ " , (66) 



2vri(t - s) Jp ^ V * - s 

and the change of variable z — (w — s)/(w — t) leads to 

T?Jx,y) = — I dw ( w-t | n e */i»-t)-v/(i»-') m (67) 

We can choose the integration path with w large and z small so that 
Re(l/(z — s) — l/(w — t)) < (in particular, z is contained in T s , so that we 
write it explicitly as r S)Z ). Then, we can exchange the integral over y with 
the integral over z and w, 

(27ri) 2 J ro 2:-+^ jT r „ (w-s)^ 1 L 



1 Jf Jf 



m-1 



(2vri) 2 J Vq z n+k+i J v ^ (w-s) n z-w' 

(6i 
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Now we enlarge the path of z so that encloses the path of w. This can be 
made at the expense of the residue at z — w. Thus we get 



(2vri) 2 J T (w-s) n J T z n+k+i z _ w 



— I dw (^~y 1 eX/{v3 - t) = ynt 



(69) 

because the first term is zero, since the residue of z at infinity is zero (k > 0). 

For the third identity, we use the representation (166]) in which we take 
the path T for w to satisfy \w\ > 1. Then, 

^•^^-^i^^T^-^ (70) 

2-ki J To (w - l)w n V / 

The last term (the integrand independent of x) is zero, because the integrand 
has residue zero at infinity, whenever n — 1 > 0. Thus, 

27T1 J ro [W — l)W n 

2ni(t - s) / ro u^ 1 J y 

= (7£ *</>)(*,</), 

The final identity is true because T™ it is the transition density of a 2n + l 
dimensional Bessel process. □ 

We proceed as in the proof of Theorem 11.21 to show that the matrix M is 
upper triangular. Indeed, with T n := l^T n t n and Lemma [3.41 we find 

c(n) ' 1 

A/r (a ,t, p ~ m W k-i\ f a n, p ~ k > tk «V{ \ J ' iii<k ' 

M k)t = {<j>** k _ t {> )(x k ) = J dxV k _ e {> (x) = < iU = k ( 72 ) 



because of the orthogonality between , k > 1, and the constant function 



Lemma 3.5. Define, for I — 0, . . . , n — 1, the polynomial o/ degree I by 

= d«, 7 -^— e-'/C-*). (73) 



2vri J r - t)«+i 
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These polynomials satisfy the orthogonal relation 

[ dx^\x)^\x) = 5 k/ (74) 
Jr + 

for k, £ = 0, . . . , n — 1. 

Proof. By the integral representation in Appendix [C] for Laguerre polynomi- 
als, we have 



t t mx/t) = — i dw- — - 

2m J Ti (w - 1 



n+£ 



and, after change of variable w = w/(w — t), we get t e Lg(x/t) = (sc) 
as defined in (1751) . The orthogonality relation (j7i]l holds because after the 
change of variable x — > xt, the left-hand side becomes 

r u\+t r 

t dxVf(xt)<f>f(xt) = • / dxxV-*L?(x)Ljf(x) = * M , (76) 
</r + \n-\-K)ii j R+ 

which is the orthogonal relation ( 1113)) for Laguerre polynomials. □ 

We now compute the kernel and start with the sum in (llOSp . Let us use 
the notations t\ = t^, t 2 = t™ 2 . Then we get 



fc=i 



If f e x x l(z-t x ) ( z _f)p-n 1 -l p "2 k 



(27ri) 2 Jr / rot2 e x ^ w -^ (w - t 2 ) p_n2+1 ~ 



We choose T and r 0i t 2 such that they do not intersect, i.e., \z\ < \w\. For 
k > n,2 the pole at w = oo vanishes and we can thus extend the summation 
over A; to oo with the result 

—L-J* dzi dw 6 Tl ^ -7 — — !— , (78) 

(2tti) 2 J ro jfr e 1 ' 2 /^-* 2 ) (u; - t 2 ) p ~ n2+1 zP w-z 



which is the second term in the kernel in Theorem 11.31 

It remains to compute ( - ta i' ta 2^. To simplify the notations, we set 
<p( ta i > ta l\x, y) = <j)( ni > tr ' n2 > t2 \x,y). By Lemma [53] we have 



V[~tT * 0* (n2_ni) , if (m.tO^Cna.ta), (7g) 



0, otherwise. 
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The integral representation fl66|) for T and (p* n (x,y) = ^t^t\ 4>(x, y) lead 
to 

^ 1 ' ii;n2 ' t2) (x ) !/) 

2vri J ro i ^ (n 2 - m - 1)! 

= Iti — t2j / dw e 

27Ti J ro M;P +1 - n 2(y; - l)T»a-ni v ' 

Finally, the change of variable u; = (z — t<i)j{z — t\) gives the first term in 
Theorem 11.31 



4 Markov property on space-like paths 

The process on matrices is clearly a Markov process along space-like paths. 
What we have to see is that the Markov property still holds for the eigenval- 
ues. The key ingredients are that the measure on matrices is invariant under 
choice of basis, and that the choice of basis at an observation point (n, t) 
depends neither on the eigenvalues at that the previous point ((n + l,t) or 
(n,t') with i! < t) nor on the eigenvalues at (n,t). 

Diffusion on GUE minors 

We first consider Dyson's diffusion. Here we denote by H(n, t) the n x n 
minor at time t and by A(n,t) the diagonalized matrix of H(n,t) which is 
obtained from conjugation by the unitary matrix U(n,t), 

H(n, t) = U(n, t)A(n, t)U*{n, t). (81) 

The Jacobian of the transformation H(n,t) — > (A(n,t),U(n,t)) gives 

dH(n, t) = A(A(n, t)) 2 dA(n, t) dfx n (U(n, *)) (82) 

where d/x n denotes the Haar measure on the unitary group U(n). 

Consider a measure at (n, t) which is invariant under unitary transforma- 
tions, i.e. with respect to the group U(n). It has the form 

A(A(n, t)) dA(n, t) d// n (£/(n, t)) (83) 

for some explicit function fi (e.g. fi(H) = exp(— Tr(H 2 )/t)). 
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Next fix n and consider times t' > t. Then, the probability measure on 
the matrices has the form (see ([35 p ) 

/i(A(n, t)) e - bTr (^(™'*')-^(™'*))) 2 dA(n, t) d^ n (U(n, t)) dH(n, t') 

= /i(A(n, t))e- bTr((A( '"'*' ) - f>( "'* )A( "'* ) ^* (ri ' t)))2 dA(n, t) d/^([7(n, t)) dff(n, 

(84) 

because of the unitary invariance of the Haar measure (here we have set 
U(n,t) = U(n,t')*U(n,t)). The integration with respect to dfi n (U(n,t)) is 
made by the well-known Harish-Chandra/Itzykson-Zuber (1115)) and the re- 
sult is as in Lemma 12.11 We are left with a probability density that depends 
only on the eigenvalues times dH(n,t'), that is, the projection onto eigen- 
values at time t did not restrict the complete freedom of choice of basis at 
(n,t'). Otherwise stated, by the decomposition (1821) . after integration over 
d[i n (U(n,t)) we have a measure on eigenvalues times dfx n (U(n,tf)) of the 
form (for some explicit f'2, which can be easily computed) 

/ 2 (A(n, t), A(n, t')) dA{n, t) dA(n, t') d[i n {U(n, t')). (85) 

The other choice is to consider t fixed and look at the measure at (n, t) 
and (n — l,t). The result explained in Proposition 4.2 of [6] is actually 
a conditional measure on eigenvalues given the eigenvalues of the minor of 
size n, thus it is not restricted to GUE, but it holds for any measure which 
is invariant under U(n), see Theorem 3.4 of [T3]. The projection on the 
eigenvalues at (n,t) and (n — l,t) leads to Lemma [2.21 We can also decide 
to project on the eigenvalues at (n,t) and (n — l,t) and the eigenvectors at 
(n— 1, t). This means that we do not integrate out the variables corresponding 
to the unitary transformations of the (n — 1) x (n — 1) minor given by 

with U(n — l,t) distributed as d/i n _i, (86) 

which form a subgroup of U(n). The eigenvalues A(n — 1, t) are independent 
of the eigenvectors (thus of the choice of basis U(n—1, t)) and the measure on 
U(n — l,t) is d/x n _i(£/"(n — l,i)) (see e.g. Corollary 2.5.4 in [3]). The measure 
on H(n, t) is invariant under U(n), so are the eigenvalues A(n, t) independent 
of the choice of U(n — l,t) (this last property follows also from the direct 
computation in Section 3.1 of [22]; Section 3.2 for Wishart matrices). Thus, 
the projection on the eigenvalues at (n, t) and (n — 1, t) and the eigenvectors 
at {n — 1, t) leads to a measure of the form 

/ 3 (A(n, t), A(n -l,t)) dA(n, t) dA(n - 1, t) d^ n _i(C/(ri - 1, t)). (87) 

for some explicit function ^3 (compare with Lemma [2.2p . 



U(n-l,t) 
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To resume, 085p and ( 187)) tell us that starting from a measure of the form 
(|83p . in which the choice of basis is completely free, the projection onto the 
eigenvalues obtained by integration over the angular variables does not fix 
the basis at the next step in the basic steps of space-like paths. This implies 
that the eigenvalues' process along space-like paths is a Markov process. 

Diffusion on Wishart minors 

Consider now diffusion on Wishart matrices. Let H(n,t) = A*(n,t)A(n,t) 
be the n x n minor at time t, where A(n, t) is the px n matrix with singular 
value decomposition A(n, t) = U (p, i)£(n, t)V*(n, t), where U(p, t) is a p x p 
Haar-distributed on U(p), V(n,t) is a n x n Haar-distributed on U(n), and 
£(n, i) is a p x n matrix with entries zeros except on the diagonal, where 
we find the singular values of A(n,t). Also, let A(n, t) = £*(n, t)£(n, t) the 
matrix of the eigenvalues of H(n,t). Thus we have 

H(n, t) = A*{n, t)A(n, t) = V{n, t)A(n, t)V*(n, t). (88) 

The Jacobian of the transformation A(n,t) — > (S(n, t), U(n,t), U(p,t)) gives 
(see e.g. [31]) 

dA(n,t) = const x (det(£*(n, t)E(n, t))) p -" +1/2 A 2 (S*(n, t)S(n, t)) 

x dS(n, t) dfi n {V{n, t)) dfi p {U{p, t)), (89) 

or, using that A(n, t) = E*(n, i)£(n, i), 

(L4(n,t) = const x (det(A(n, t))) p ~ n A 2 (A(n, t)) 

x dA(n, t) d/jL n (V(n, t)) d/jL p (U(p, t)). (90) 

Therefore, the starting measure at (n, t) has the form 

gx(A(n, t)) dA(n, t) dfi n (V(n, t)) dfi p (U (p, i)) (91) 

for some explicit function g\. 

Next consider fixed n and time t' > t. Then, the probability measure on 
the matrices has the form (see (1571) ) 

0l(A(n, *)) e - 6Tr ((^*(n/)-^*(M))(^(™/)-^(n,t))) 

x dA(n, t) d/i n (\/(n, t)) dn P (U(p, t)) dA(n, t') 

= gi(A(n, t)) e -6Tl'(P*(n.*')-V*(n»t)S , (n,t)&(p,t)] [E(n,f )-^(p,t)E(n,t)V-(n,t)]) ( 92 ) 

x dA(n, t) d/i n (\/(n, t)) dp, p (U(p, t)) dA(n, if) 
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because of unitary invariance of the Haar measure (we set V(n, t) = 
V(n,t')*V(n,t) and U(p,t) = U(p,t')*U(p,t)). An integration with re- 
spect to dfi n {V(n, t)) d/i p ([/(p, t)) according to (I116P results in the formula 
of Lemma 13.11 We are left with a probability density that depends only 
on the eigenvalues times dA(n, t'), that is, the projection onto eigenval- 
ues at time t did not restrict the complete freedom of choice of basis at 
(n, t'). Otherwise stated, by (MI]) we have a measure on eigenvalues times 
d/i„(V(ri, t')) dn P {U{p, t')) of the form (for some explicit which can be 
easily computed) 

g 2 (A(n, t),A(n, t')) dA(n, t) dA(n, t') dfi n (V{n, t')) d/i p (U(p, t')). (93) 

The other choice is to consider t fixed and look at the measure at (n, t) 
and (n — l,t). This works as for the Hermitian case and we get a measure 
of the form 

g 3 (A(n, t), A(n - 1, t)) dA(n, t) dA(n - 1, t) d^iV^n - 1, *)) d/i p (U(p, t)). 

(94) 

for some explicit function g% (compare with Lemma [3. 2ft . 

Therefore (1931) and (194^) tell us that starting from a measure of the form 
(I9ip . in which the choice of basis (in which the matrix A is represented) is 
completely free, the projection onto the eigenvalues obtained by integration 
over the angular variables does not fix the basis at the next step in the basic 
steps of space-like paths. This implies that the eigenvalues' process along 
space-like paths is a Markov process. 

A Space-like determinant al correlations 

For convenience we report here the statement of Theorem 4.2 of [9]. 

Let Xi, . . . , £jv be finite sets and c(l), . . . , c(N) be arbitrary nonnegative 
integers. Consider the set 

X = (Xi U • ■ ■ U Xi) U • • ■ U (X N U ■ ■ • U X N ) (95) 

with c(n) + 1 copies of each X n . We want to consider a special form of the 
weight W(X) for any subset X C X, which turns out to have determinantal 
correlations. 

To define the weight we need a bit of notations. Let 

4> n ( ■ , • ) : X n _i xI„4C, n = 2,...,N, 

n (virt, • ) : X n ->> C, n = l,...,N, (96) 
v&f(-):Xjv^C, j = 0,...,iV-l, 
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be arbitrary functions on the corresponding sets. Here the symbol virt stands 
for a "virtual" variable, which is convenient to introduce for notational pur- 
poses. In applications virt can sometimes be replaced by +00 or —00. The 
(fi n represents the transitions from £ n _i to X n . 
Also, let 

£() — — l c(N) — l ^ — l c(N-l) ~ l — - r c(2) — r ^ ^ l c(l) 

(97) 

be real numbers. In applications, these numbers refer to time moments. 
Finally, let 

7*S>*"-i( ' ' ' ) : %n x £„ -> C, n = l,...,N, a = 1, . . . ,c(n), (98) 

be arbitrary functions. The 7t»,t« represents the transition between two 
copies of X n associated to "times" t™_ x and t™. 

Then, to any subset Iclwe assign its weight W(X) as follows. W(X) 
is zero unless X has exactly n points in each copy of X n , n = 1, . . . , N. In 
the latter case, denote the points of X in the mth copy of £„, by x£(t^), 
k — 1, . . . , n, m — 0, . . . , c(n). Thus, 

X = K(C) I fc = 1, • • • , n; m = 0, . . . , c(n); n = 1, . . . , N}. (99) 

Set 



iV 



h/(x) = n 

n=l 

c(n) 

X 



det [^(^r 1 ^- 1 ),^^)))] 



(n)^J Kfc ;<n 



Hdet [^(xjco.xrcci))]^,, 

a=l 



where x ( • ) = virt for all n = 1, . . . , N. 



det [K-^ntm^,^ 

(100) 



In what follows we assume that the partition function of our weights does 
not vanish, 

Z :=Y^ W ( X ) ± °- ( 101 ) 



xcx 



Under this assumption, the normalized weights W(X) = W(X)/Z define a 
(generally speaking, complex valued) measure on 2 X of total mass 1. One 
can say that we have a (complex valued) random point process on X, and its 
correlation functions are defined accordingly, see e.g., [13]. We are interested 
in computing these correlation functions. 
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Let us introduce the compact notation for the convolution of several tran- 
sitions. For any n — 1, . . . , N and two time moments t™ > t% we define 

Tt n in = I in m * lf n t n * ■ ■ ■ * Tf n )», / n = Tt n f n , (102) 

where we use the notation (/ * g)(x,y) := J2 Z f( x > z )di z i v)- For any time 
moments t™* > with (ai,ni) 7^ (02,^2), we denote the convolution over 
all the transitions between them by (jp^^i; 

(t ai ,* a2 ) = T + + T n 1+ 1 + . . . ^ 7- (1Q3) 

If there are no such transitions, i.e., if i™* < or (ai, n\) = (02, rc.2), we set 
Furthermore, define the matrix M = || Mk t i\\^ l=1 by 

M M = (fa * T fe * ■ • • * <j> N * T N * K-i) (virt) (104) 

and the vector 

<5 = ^ * / = l,...,iV. (105) 

The following statement describing the correlation kernel is a part of Theorem 
4.2 of 0. 

Theorem A.l (Part of Theorem 4.2 of [9]). Assume that the matrix M 
is invertible. Then Z = det M 7^ 0, and the (complex valued) random point 
process on X defined by the weights W(X) is determinantal. Its correlation 
kernel can be written in the form 



N n 2 

+ EEC 1 W[ jr Wi * { ^ A2) )(virt,x 2 ). (106) 



k=l 1=1 



Remark A. 2. As stated in the complete statement of Theorem 4.2 of [9], 
there is one situation where the kernel takes a simple formula. Namely, when 
the matrix M is upper-triangular, then 

®2-H x ) ■= Y). M ~\i(<l>i * ( ' iwA2) )(virt,x) (107) 



1=1 



are the function biorthogonal to ^ n ^k (#) obtained for the non-extended 
kernel (i.e., at fixed level and fixed time). In the case of random matrices 
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which we consider, the functions , k = 0, ... ,n — 1, have to be poly- 
nomials of degree k because det {xj))i<j,k<n must be proportional to 
A(x) = det(x^~ 1 )i<j i fe< n , the Vandermonde determinant. Then, the kernel is 
simply written as 

K{n l X\^ l] n 2 Xlx 2 ) = -^^ 

k=l 

(108) 

Remark A. 3. Looking at the proof of the above theorem in [9] one also sees 
that the time evolutions T can be taken to be level-inhomogeneous, i.e., the 
77« a n can be a function of i",in and also of the level n. Such a situation 
occurs for the Wishart matrices case. 

The proof of Theorem IA.1I given in [S] is based on the algebraic formalism 
of [13]. Another proof can be found in Section 4.4 of [21]. Although we stated 
Theorem IA.1I for the case when all sets X n are finite, one easily extends it to 
a more general setting. Indeed, the determinantal formula for the correlation 
functions is an algebraic identity, and the limit transition to the case when the 
X n 's are allowed to be countably infinite is immediate, under the assumption 
that all the sums needed to define the *-operations above are absolutely 
convergent. 



B Hermite polynomials 

The Hermite polynomial of degree n is denoted here p n (x). We use the 
normalization of [31] : 

/ dxe-' x2 p n (x)p m (x) = b m ^Tn\. (109) 
Jr 

There are two useful integral representations for the Hermite polynomials 



on 

1 ^ 



p n [x) = -^e x / dwe w - zxw w n : 
Pn {x) = — l dze- (z2 - 2xz h-( n+1 \ 



2vri 



;no) 



(i 
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as well as the identities (with < q < 1) which can be found in 



(x - qy) 2 \ _ x 2 Pk(x)Pk(y)q 



exp (-^r^- =e £ 



-,k 



in: 



2 2 

dye" 5 ' p n (j/) = e - * p„_i(a;), 

p n (x) = (-l)> n (-x). 

These identities can be useful to rewrite the double integral representation 
into an expression in terms of Hermite polynomials (as it was made e.g. in 
Lemma 24 of [12] for the antisymmetric GUE minor kernel). 



C Laguerre polynomials 

The generalized Laguerre polynomials L\ of degree k and order p are poly- 
nomials on M + defined by 

They satisfy the orthogonal relation 

/ dxx p e- x Ll(x)L p e (x) = ^-^-5 k/ (113) 
Jr+ fc! 

and have integral representations, 



i r p -x{w-l) p+k 



2vri fr, (w - l) k+1 



(114) 



klxP 2ni J To z p+ k+1 



D Harish-Chandra/Itzykson-Zuber formulas 

Here we report the Harish-Chandra/Itzykson-Zuber formula as well as its 
generalization for rectangular matrices. 

Let A = diag(ai, . . . , ajv) and B = diag(&i, . . . , two diagonal N x N 
matrices. Let d/i denote the Haar measure on the unitary group U{N). Then, 



det (e^) K . . <N 
dfi(U)exp(Tr(AUBU*)) = J[p\ 1(a) A(6) ' ^ 
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where A (a) is the Vandermonde determinant of the vector a = (a 1; . . . , a^). 

The extension to rectangular matrices can be found in section 3.2 of (12] 
and was derived in [25]. Let A be a complex N\ x matrix, B a complex 
./V2 x N\ matrix so that the A2 x A2 matrices A* A and BB* are diagonal with 
(real positive) eigenvalues a = (ai, . . . , a^ 2 ) and b = (b±, . . . , 6jv 2 ) respectively. 
W.l.o.g. we assume iV\ > N 2 . Then, 



dn(U) / dfi(V)exp(Tr(AUBV* + B*U*A*V)) 

U{N 2 ) JU(Ni) 

_ ngrV nr=r x j det yj^V^^^ 

n ^-Ar 2 -i H A(fl)A(6) n£(a,^)W-^)/ 2 ' 
where 7 n is the modified Bessel function defined by 



(116) 



1 r p x(z+z x ) 00 ™fc ™fc+|n| 
IJ2X) = — (b dz- -— = Yt T /1 , [TT ) (117) 

2vri7 ro z^ 1 ^ jfe! (jfe + |n|)!' V ; 



for n G Z. 
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